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Abstract 

We have developed the method that allows us to estimate 
the magnetic field strength at the horizon of a supermassive 
black hole (SMBH) through the observed polarization of op- 
tical emission of the accreting disk surrounding SMBH. The 
known asymptotic formulae for the Stokes parameters of out- 
going radiation are azimuthal averaged, which corresponds to 
an observation of the disk as a whole. We consider two mod- 
els of the embedding SD-magnetic field, the regular field, and 
the regular field with an additional chaotic (turbulent) com- 
ponent. It is shown that the second model is preferable for 
estimating the magnetic field in NGC 4258. For estimations 
we used the standard accretion disk model assuming that the 
same power-law dependence of the magnetic field follows from 
the range of the optical emission down to the horizon. The 
observed optical polarization from NGC 4258 allowed us to 
find the values 10'^ — 10* Gauss at the horizon, depending 
on the particular choice of the model parameters. We also 
discuss the wavelength dependencies of the light polarization, 
and possibly applying them for a more realistic choice of ac- 
cretion disk parameters. 

Keywords: polarization - magnetic fields - accretion disks 
- supermassive black holes; galaxies: active. 

1 Introduction 

It is now commonly accepted that active galactic nuclei 
(AGNs) and quasars (QSOs) frequently possess the magne- 
tized accretion disks (see, for example, the reviews of Blaes 
2003; Moran 2008 on NGC 4258). There are many models of 
the accretion disk structures (see Pariew et al. 2003, and ref- 
erences therein) . The best known and most frequently used is 
the standard model of Shakura & Sunyaev (1973). The polari- 
metric observations frequently demonstrate that AGNs and 
QSOs have polarized emission in different wavelength ranges, 
from ultraviolet to radio waves, in continuum and in the line 
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emission (see Martin et al. 1983; Webb et al. 1993; Impey 
et al. 1995; Wilkes et al. 1995; Barth et al. 1999; Smith 
et al. 2002; Modjaz et al. 2005). These papers discuss the 
different mechanisms for the origin of the observed polariza- 
tion: the light scattering in accretion disks, which happens 
on both free and bound electrons, synchrotron radiation of 
charged particles. These mechanisms can work in different 
structures such as the plane and warped accretion disks and 
toroidal clumpy rings, surrounding the accretion disks and 
jets. Frequently different models are proposed to explain the 
same source. There are a lot of papers devoted to different 
aspects of the structure and emission of AGNs and QSOs. 
Many theoretical papers propose the possible behavior of a 
magnetic field in these objects. 

In this paper we develop the technique of estimating the 
magnetic fields in different parts of plasma accretion disks. 
Especially interesting is the estimation of magnetic field in 
the horizon of the supermassive black holes in AGNs. The 
main idea is to use the observed integral polarization from 
magnetized plasma accretion disks. We use the known fact 
that the Faraday rotation of polarization plane changes both 
the values of integral polarization degree p and position angle 
X- The observed spectra p{\) and x(A) acquire very specific 
forms due to Faraday rotation. The detailed discussion and 
calculations of these effects are presented in Silant'ev (1994), 
Dolginov et al. (1995), Gnedin & Silant'ev (1997), Agol and 
Blaes (1996), etc. The observed polarization possesses the 
information about the magnetic field in magnetized electron 
atmospheres and can serve for estimating the field. In Gnedin 
et al. (2006) the method was considered for pure vertical mag- 
netic field B|| and without the correct azimuthal averaging of 
the asymptotic formulae (4). 

For estimating the magnetic field we use the simple ap- 
proximate formulae (Silant'ev 2002) that represent solutions 
to a number of "standard" problems of the radiative trans- 
fer theory in magnetized electron atmospheres, namely, the 
Milne problem and the cases when the sources of thermal 
radiation are distributed homogeneously, linearly, and expo- 
nentially in an optically thick atmosphere. These "standard" 



1 



N Z 




Figure 1: Main notions of accretion disk surface, geometry of 
magnetic field B, and line of sight n (see text). 



solutions allow us to approximate the solution of problem 
with a more complex distribution of thermal sources inside 
the atmosphere, because the latter can be presented as a su- 
perposition of "standard" sources. 

For the optically thick accretion disks, the solution of the 
Milne problem is used, i.e. the case where the sources of 
thermal radiation arc located far from the surface. The po- 
larization and angular intensity distribution of outgoing radi- 
ation for non-magnetized electron atmosphere is presented in 
the known Chandrasckhar's book (sec Chandrasckhar 1950). 
The numerical solution to Milne's problem for a magnetized 
electron atmosphere with the magnetic field B|| parallel to 
the normal N to an atmosphere is presented in Agol & Blacs 
(1996), and Shternin et al. (2003). The numerical solution 
to this problem for the turbulent magnetized atmosphere is 
given in Silant'ev (2007). In this paper the approximate for- 
mulae for Milne's problem are also generalized for the case of a 
turbulent atmosphere. It is very important that the approxi- 
mate formulae of Silant'ev (2002, 2007) are valid for arbitrary 
directed magnetic field. 

Some words would be useful here about the simplifications 
in our method. First of all, we consider the optically thick 
plane plasma accretion disk neglecting the possible warps. 
Using the Milne problem we neglect the reflection of radi- 
ation from possibly existing central outflows (frequently this 
radiation lies far from observed optical wavelength bands). In 
all models of magnetized accretion disks, the solutions within 
the framework of the power-law dependence of magnetic fleld 
are sought inside the disk (see, for example, Pariev et al. 
2003). Physically this assumption seems fairly natural if we 
remember that far from the sources the magnetic fields tend 
to dipole, quadrupole etc. forms, i.e. acquire the power-law 
dependence. We also assume that the radial dependence of 
the disk's magnetic field follows the same power-law in the 
range from the optical polarized emission down to the hori- 
zon. These simplifications now are commonly accepted, and 
can be considered as important assumptions of our theory. 



The Milne problem in terms of vertical Thomson depth 
includes the possible vertical inhomogeneities of the atmo- 
sphere, so we do not include only possible horizontal inho- 
mogeneities of the atmosphere. But if these inhomogeneities 
are smooth (with the characteristic length of many Thomson 
free lengths), the corrections should be neglected. In our pa- 
per we do not include the true absorption effects, considering 
the Milne problem in the limit of conservative atmosphere. 
Certainly, some our simplifications, such as the latter one, 
can easily be taken into account in the proposed method. 
We stress that this method can be generalized to more com- 
plex situations; in particular, it may be considered together 
with the other sources of polarized radiation (polar outflows, 
toroidal clumpy disks, etc.). 

It should be mentioned that many AGNs models postulate 
the existence of a dusty geometrically thick obscuring region 
"the torus," which is placed far from the center of AGN (see 
Chang et al. 2007, and many references therein). This region 
give additional infrared radiation, as compared to the usual 
radiation of the interstellar medium. The spectrum of linear 
polarization from dusty media is characterized by Serkowski's 
formula (see Serkowski 1973; Martin 1989). If the spectrum 
of polarization of an AGN differs strongly (as in the source 
NGC 4258) from Serkowski's distribution, then the probabil- 
ity that the polarization comes from multiple scattering in 
plasma disk increases. 

Our goal in this paper is to present the method of estimat- 
ing of magnetic fields for fairly simple models. For this reason, 
in particular calculations we restrict ourselves to the most 
popular standard disk model of Shakura & Sunyaev (1973). 

2 Basic equations 

We begin with the known expression for the Faraday rota- 
tion angle at the Thomson optical path r = arZ that is 
frequently used below: 
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where n± are the refractive indices for right and left circular 
polarized electromagnetic waves, A = 2-kc/w the wavelength 
of the radiation and 6 the angle between the magnetic field B 
and line of sight n directions, or Thomson cross-section, and 
Ng is the number density of electrons. The first expression in 
Eq. (1) is presented in many textbooks on optics. The second 
expression is derived in books on plasma physics (e.g. KroU 
& Trivelpiece 1973). This expression is usually used in radio 
astronomy (see Rohlfs & Wilson 1996). In radiative transfer 
of optical radiation, the third expression in Eq.(l) is prefer- 
able, because it uses well known Thomson optical depth r (see 
Dolginov et al. 1995, Gnedin & Silant'ev 1997). According to 
Eq. (1) the Faraday dimensionless depolarization parameter 
5 takes a simple form: 
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Here to = 2'kv = 27rc/A is angular frequency, ujb = eB/nieC 
the cyclotron frequency of an electron in a magnetic field, 
Te = e'^/niec'^ — 2.82 • 10~^^ cm is the classic electron radius. 
Parameter ujb/<jJ — 0.93 • 10^^ X{fim)B{G) is assumed to be 
small. This is the condition where the simple Eq. (1) is valid. 
In our paper we have B « 10^ — lO^G near the black hole's 
horizon. It this case the possible wavelength A can acquire 
values from X-rays up to infrared radiation with A ~ lOO/xm. 
Below we use the asymptotic formulae for the Stokes param- 
eters, which arc valid for > 1. Therefore, the wavelength 
range for which our model is valid extends from the X-ray 
band up to A ~ 0.1/um. Far from the horizon the magnetic 
field is less, and the depolarization parameter 5-^1. In this 
case we have to use the usual Chandrasekhar formulae for 
polarization. 

Silant'ev (2002) derived the asymptotical analytical formu- 
lae for the Stokes parameters of the radiation emitted from 
a magnetized, optically thick, plane-parallel atmosphere. For 
the Milne problem they are 
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where /7, = cos I is the cosine of the angle between the normal 
to atmosphere N and the line of sight n, q the degree of 
true absorption {q = Uaj {pa + fTs)),and F\ the total radiation 
flux. Function J(/i) describes the angular distribution of the 
radiation emerging from a disk. This function, as well as the 
numerical parameters g, fc, and J\ were tabulated by Silant'ev 
(2002). For electron conservative atmosphere (g — 0), the 
values of these parameters are A: = 0, 3 = 0.83255, and J\ = 
1.19402. The Stokes parameters Q and U are given in the 
reference frame with the X-axis lying on the plane (nN) (see 
Fig. 1). 

Formulae (3-4) for polarization consider the last scattering 
of radiation before escaping from a semi-infinite magnetized 
atmosphere. For a high value of parameter (5, the contribu- 
tion of the secondary scattered photons is small because of 
the large Faraday depolarization. Even in the absence of a 
magnetic field, the main contribution to the polarization of 
emitted radiation comes from the last scattered photons. For 
this reason, Eqs. (3-4) at the absence of magnetic field practi- 
cally represent the classical Chandrasekhar- Sobolev polariza- 
tion in the Milne problem (see, for example, Chandrasekhar 
1950). 

As far as the intensity of radiation /(/x) is concerned, one 
can remember (see Chandrasekhar 1950) that the polarization 
weakly influences the intensity. For Milne's problem without 
the true absorption (g = 0), we have /(O) ~ 3.06, whereas the 
separate transfer equation with the Rayleigh phase function 
gives 7(0) ~ 3.02. For high values of 5, the terms with Stokes 



parameters Q and U in the full system of transfer equations 
for parameters /, Q and U become very small ~ 1/5, and they 
are negligible in the equation for intensity I. As a result, 
the radiation intensity obeys the separate transfer equation 
with the Rayleigh phase function (see, Silant'ev 1994 for more 
detail). Expression (3) presents the solution to this equation. 
For large 8 the main contribution to polarization comes from 
to intensity term. Formulae (4) were obtained in this way. 

Equations (2-4) allow us to derive the following approxi- 
mate expressions for polarization degree p(B, n) and the po- 
sition angle x of radiation for an accreting magnetized disk: 
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where B cos 6 = Bn. Now we consider the case of dominant 
nonabsorbing electron scattering in accretion disks, i.e. g = 
and fc = 0. Photons escape the optically thick disk basically 
from the surface layer with r w 1. If the Faraday rotation 
angle ^ corresponding to this optical length becomes larger 
than unity, then the emerging radiation will be depolarized 
as a result of summarizing the radiation fluxes with the very 
different angles of Faraday rotation. Only for directions that 
are nearly perpendicular to the magnetic fleld are the Faraday 
rotation angles too small and depolarization does not occur. 
Certainly, the diffusion of radiation in the inner parts of a 
disk depolarizes the light, even in the absence of a magnetic 
fleld, because of multiple scattering of photons. 

The existence of a magnetic field, hence Faraday rotation, 
only increases the depolarization process. It means that the 
polarization of outgoing radiation acquires a peak-like angular 
dependence with its maximum for perpendicular propagation. 
The sharpness of the peak increases with increasing magnetic 
fleld magnitude. The main region of allowed angles appears 
to be ~ 1/5. 

Another very important feature of the polarized radiation 
is the wavelength dependence of polarization degree p and po- 
sition angle x which is very different from the case of classical 
electron scattering. This effect is briefly considered in Sect. 5. 
We consider the case where the Thomson cross-section does 
not depend on the radiation wavelength. 

2.1 Integral polarization from the accretion 
disk 

The axially symmetric accretion disks frequently are observed 
as a whole. The observed integral Stokes parameters (Q) and 
IJJ) are described by the azimuthally averaged formulae (3) 
and (4). To derive these expressions we introduce the follow- 
ing notions: 



5cos6 = 5|| cosi + 6± sini cos(<^ -|- tp*) = a + fecos$, (6) 

where ip is the azimuthal angle of radius-vector r (the az- 
imuthal angle of line of sight n is taken zero), (p» the an- 
gle between perpendicular magnetic fleld projection Bj_ = 
{Bp, By,), and the projection Bp (cos = Bp/Bj^) (see 
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Fig. 1). The value By = is the magnetic field projection 
along the normal N to the accretion disk plane. According 
to Eqs. (1) and (6) we have 

a = di\fi, = 0.8A2(/zm)S|| (G), 

b = S±^/l-^J.^, S_L = 0.8X'^{nm)B_LiG). (7) 

Using the axial symmetry of an accretion disk, we obtain (re- 
member that U (0, /x) = for Milne problem without magnetic 
field): 

{Q) = Q{0,^^)x 

l + a^+6^cos^^ 
(1 + a2 + 62 cos2 $)2 - {2ab cos $)2 ' 
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The observed degree of the light polarization and the position 
angle are derived from parameters (8) in the usual way. For 
particular cases of pure normal {6j_ = 0) and pure perpendic- 
ular = 0) magnetic fields, expression (8) can be derived 
analytically. For the first case we have 

p(B,n) = ^B^, tan2x = ^||/x- (9) 
For perpendicular magnetic field are there the formulae 
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For a pure perpendicular magnetic field, the position angle 
X = comes from the axial symmetry of the problem, so 
the electric wave oscillations in this case occur parallel to the 
surface of an accretion disk. The case x 7^ can be realized 
if some Bz component exists. It is interesting that the cases 
{Bp ^ 0,B^ = 0) and {Bp = 0, ^ 0) give rise to the 
same formula (10). This is the consequence of the averaging 
procedure. 

It is seen from Eqs. (8) that the relative degree of polar- 
ization p(B, n)/p(0, /i) and position angle x only depend on 
dimensionless parameters a and 6, which are the functions 
of wavelength A, magnetic fields B|| or B±, and inclination 
angle i (/i = cosi). First we discuss the behavior of the rel- 
ative degree of polarization and position angle on these two 
parameters. Remember that we consider conservative atmo- 
sphere with q = 0, hence k = 0. It is interesting to investigate 
how the polarization changes if we include the perpendicular 
magnetic field B± in the existing parallel magnetic field B|| . 

The average process takes the disk regions with the very dif- 
ferent angles 6 into account between the magnetic field B and 
the line of sight n, and the integral polarization and position 
angle can acqiiire very different values. The numerical calcu- 
lations show that position angle x not only depends on param- 
eter a = but also on the 6 - parameter (6 = 5j_-^l — /jp). 
The existence of perpendicular magnetic field Bj^ diminishes 



the value of x compared to the case of pure parallel magnetic 
field. This decrease is especially large if 6 > a and the pa- 
rameters a and 6 are close to unity. For a ^> 1 and b < a 
the decrease; of x is small and practically x — 45°. But for 
6 > a » 1, the position angle x 0. The special case is 
a = 6 3> 1. In this case the position angle x rapidly decreases 
from the limiting value ~ 45° for a — 6~ 1-^5 to value 22.5° 
at 6 = a, and then tends to zero for 6 — o ~ l-j-5, so for large a 
the intermediate values of x can only occur in a rather narrow 
interval 6 — a ~ 1 2, i.e. at a ~ 6. 

The degree of linear polarization p depends on parameters 
a and 6 in a more complex form. For a < 1 the addition of the 
perpendicular magnetic field (parameter 6) lowers the integral 
polarization. For a > 1 there is the region of 6 (6 < a) where 
the polarization increases compared with the case of pure par- 
allel magnetic field. The maximum polarization occurs at 
6 ~ a, and then the polarization decreases with the increase 
in parameter 6. The increase in polarization at h = a can 
be rather large. As a result, for value a = 6 = 5, 10, 20, and 
50 the relative polarization increase, as compared to purely 
parallel magnetic field, is equal to 160%, 224%, 317%, and 
504%, respectively. It seems this effect stems some "reso- 
nant" regions in an accretion disk where the Faraday rotation 
from parallel magnetic field is balanced by opposite rotation 
from a perpendicular magnetic field. Of course, the magni- 
tude of polarization decreases with the increase in a and 6. 
The numerical calculations demonstrate that the relative po- 
larization degree p{'B,n)/p{0, n) is a symmetric function of 
parameters a and 6. The position angle x does not possess 
this symmetry. 

Now we shortly discuss the wavelength dependence of po- 
larization degree p{X) and xi^)j which follows from general 
formulae (8). More detailed discussion is presented in Sect. 
5. For high values of parameters a and (or) 6 the spectra di- 
minish ^ 1/A2. But for the case a = b mentioned above, the 
spectra diminish as ^ 1/A. Thus, the "resonant" effect also 
changes the asymptotic behavior of spectra. Note once more 
that this effect disappears beyond the interval |ci — 6| « 1 — 5. 
The spectra x(-^) depend strongly on the relative value of 
the perpendicular magnetic field B± as compared to vertical 
component -B|| . If parameter a » 1 and a ^ b the position 
angle x(A) — > 45°; i.e., it becomes independent of wavelength. 
For the "resonant" case a = b, this limiting value is equal to 
22.5°. 

The characteristic spectra of polarization and position an- 
gle as a function of inclination angle i and magnetic fields 
and B± are presented in Fig. 2 (/i = cosi). The men- 
tioned case a = b corresponds to B^^fj, = B^\/l — jj?] i.e., at 
B\\ = B± it exists at i = 45°. For every A there exists its own 
value a(A) = 6(A), so, for A = 1/im and B|| = B± = 5G, this 
value is a = 6 = 2.828. In this figure we present the relative 
polarization degree p(B, n)/p(0, /u) for the inclination angles 
i = 85°, 60°, 45°, and 30°. The values of the polarization de- 
gree p{0, jj) for these angles are equal to 7.80%, 2.25%, 1.08%, 
and 0.43%, respectively (see Chandrasekhar 1950). The pre- 
sented spectra can help readers recognize general tendencies 
of polarization as a function of the basic system parameters. 

We see that the higher magnetic field, the larger depolariza- 
tion. For i K, 90° and pure parallel magnetic field (Bj^ = 0), 
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Figure 2: Spectra of relative polarization degree p(B,n)/p(0,/j) and position angle x- The numbers denote the values of inclination 
angle i. The numbers in brackets refer to dotted curves. The a) figures demonstrate spectra for the pure parallel magnetic field with 

B|| = 5 G (bold curves) and 50 G (dotted curves). The b) figures concern the case B\\ — B± = 5 G (bold curves) and 50 G (dotted 
curves). The c) figures present spectra for B|[ = 20 G and B± = 10 G (bold curves), and for By = 10 G and B± = 20 G (dotted 
curves) . 

the polarization degree tends to the Thomson value of polar- the magnetic field is practically perpendicular to the line of 
ization p(0, fi). This is quite natural because in these cases sight n, and the Faraday rotation is low. The position angle x 
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is more sensitive to Faraday rotation, and tends to Thomson 
value (x = 0) slower than the polarization degree tends to the 
Thomson polarization. It is interesting that for this case the 
relative polarization degrees p(B, n)/p(0, /x) and position an- 
gles X practically coincide for inclination angles i = 85° and 
i = 30°, if the magnetic fields differ 10 times (for example, 
S|| = 5 and 50, or B\\ = 10 and 100, etc.). This happens 
because the corresponding values of = cosi, 0.08715 and 
0.86602, differ approximately 10 times. 

Comparing formulae (9) with (10), we find that the relative 
polarization degrees for perpendicular magnetic field (i?|| = 0) 
can be taken from the results drawn in Fig. 2a, if one uses 
the substitution i (90° — i) there. The case B± = 5 G 
and i = 60° therefore coincides with the case = 5 G and 
i = 30°. Of course, the position angle x = for perpendicular 
magnetic field. 

For the case = B± (see Fig. 2b) the most depolarization 
occurs at z ~ 90°, in contrast to the pure parallel magnetic 
field. It is interesting that the relative degrees of polarization 
are the same for the inclinations i and 90°— «. But the position 
angles are different in these cases. The relative polarization 
degree is higher for i = 45° than for i = 30°, i.e. the change 
in this value is not monotonic. The case i = 45° corresponds 
to equality a = b ^ 1. As mentioned above, in this case 
X — >■ 22.5°. The righthand side of Fig. 2b confirms this. 

Figure 2c presents the spectra of the relative polarization 
degree and position angle for B|| = 20 G and B± = 10 G (the 
curves denoted by the usual numbers) and the opposite case 
= 10 G and B± = 20 G (the curves denoted by numbers 
in brackets). The symmetry of relative polarization degree as 
a function of parameters a and b gives rise to the coincidence 
of spectra in the first case with those in the second case if the 
angle i — > (90° — i), so the relative polarization degree spectra 
at 30°, 45°, 60° and 85° coincide with the spectra denoted as 
(60°), (45°), (30°) and (5°), respectively. But the position- 
angle spectra are diflferent for these two cases. 

The spectra, presented in Fig. 2, demonstrate a large vari- 
ability of values and forms for polarization degrees and posi- 
tion angles in the integral radiation escaping from the mag- 
netized accretion disks. 

2.2 Some results from accretion disk models 

Models of a magnetic accretion disk with externally imposed, 
large- scale vertical magnetic field and anomalous magnetic 
field diffusion due to enhanced turbulent diffusion have been 
considered by Campbeh (2000), Ogilvie & Livio (2001), and 
Pariev et al. (2003). 

We calculate the value of Faraday depolarization parameter 
6 for the model of an accreting disk suggesting the power-law 
radial dependence of the magnetic field: 

B{r)=BH{RH/rr, (11) 
where B{r) is the magnetic field inside an accretion disk, and 
Bh the magnetic field strength at the event horizon radius of 
a SMBH. The radius of the black hole horizon is (Novikov & 
Thorne 1973) 

Rh = (l + Vl - (a/MBHr) , (12) 



where a/MsH is the spin of a Kerr black hole. 

Pariev et al. (2003) have developed the detailed descrip- 
tion of magnetized accretion disks with the different values for 
parameter n (denoted as S in their paper). The radial effec- 
tive temperature dependence Te ~ r~^/^ in their model is the 
same as in the Shakura-Sunyaev model. They find that, in the 
case of equipartition of magnetic pressure with radiation or 
thermal pressures, their results are close to Shakura-Sunyaev 
model with the viscosity parameter a = 1. As a most phys- 
ically significant model, they investigate the case n = 5/4 in 
more detail. Parameter n must be greater than unity if we 
consider the disks with the diminishing gas density far from 
the black hole. It should also be noted that dependence (11) 
with n = 5/4 takes place at the equipartition of magnetic 
energy with thermal and gravitational ones in the spherical 
accretion (see, for example, Melia 1992). Following to Pariev 
et al. (2003), we use mostly the case n = 5/4. 

The central problem is to derive the characteristic scale R\ 
that corresponds to the effective wavelength of polarimetric 
observations. At first glance, we can estimate the radius Rx, 
suggesting that A = Ae corresponds to the value of rest-frame 
wavelength of the black body spectrum maximum, 

29 

Ae = -— , Ao = Ae(l + z), (13) 
Te{r) 

where Aq is the wavelength in the observer system, and z the 
cosmological redshift. 

In a standard thin disk model (Shakura & Sunyaev 1973), 
there are a black body radiation with an effective temperature 
profile of Te = TH{r/ Rh)~^^^ and the scale length R\ defined 
by the point in the disk where the disk temperature matches 
the rest-frame wavelength of the monitoring band. 

There is the series of papers where the semi-empirical 
method of determining of the accretion disk scale R\ has been 
developed (see Kochanek et al. 2006; Poindexter et al. 2008; 
Morgan et al. 2007, 2008). The authors used microlensing 
variability observed for gravitationally lensed quasars to find 
the accretion disk size and the observed (or rest-frame) wave- 
length relation. It is very important that the scaling appeared 
to be consistent with what is expected from the thin accretion 
disk theory of Shakura & Sunyaev (1973). 

This allows us to have the following size scaling (Poindexter 
et al. 2008): 

«^(»')— °-(^)"(^)"(?)"' (») 

Here R\ is the distance in the accretion disk, which corre- 
sponds to Ae; ?7 = Lhoi/LEdd, Lsdci = 1-3 x 10^^{Mbh/Mq) 
erg s~^ as the Eddington luminosity, and e as the rest-mass 
radiation conversion efficiency. Note that R\ does not depend 
on viscosity parameter a. We use the commonly accepted re- 
lation between the bolometric luminosity Lboi of the accretion 
disk and the accretion rate M: 

Lboi = eMc^. (15) 

Theoretical calculations of parameter e depend on the de- 
tails of the accepted model. Usually one presents the depen- 
dence of e on the spin parameter a/MBH for every particular 
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model. Naturally, this dependence is different for different 
models. Krolik (2007) present the comparison of £-values 
from Novikov-Thorne model (Novikov&Thorne 1973) and the 
model that take the jet luminosity into account. For spin 
parameter a/MBH = 0.5, the first model gives e = 0.081, 
whereas the second model gives e = 0.0063. In some papers 
one considers s as an independent parameter in the model 
calculations. For example, Pariev & Colgate (2007) accept 
£ = 0.1, rj = 0.1 and viscosity parameter a = 0.01. 

The strong gravitational field near the black hole infiuences 
the Stokes parameters of outgoing radiation when they prop- 
agate to an observer. The detailed calculations of this effect 
have been done by Connors ct al. (1980), Karas et al. (2004), 
and Dovciak et al. (2004). Usually these effects are impor- 
tant when describing X-ray emission from the vicinity of the 
black hole. The optical radiation arises far from this place, 
and we neglect these gravitational corrections. 

3 Magnetic field strength of NGC 
4258 

NGC 4258 is a low-luminosity Scyfcrt II galaxy at the distance 
of about 7.2 Mpc, which harbors water masers (Modjaz et 
al. 2005). This object is usually considered as a very good 
laboratory for successfully measuring the magnetic field in 
accretion disk even very close to the central black hole. The 
spatial velocity distribution of water mega-maser sources in 
NGC 4258 on scales of 0.14-0.28 pc indicates a thin Keplerian 
disk rotating around a black hole with a mass Mbh = 3.9 x 
1O'^M0 (Herrnstein et al. 1999). The accretion disk has an 
almost edge-on orientation with the radiation axis and its 
inclination angle is « = 83° ± 4° (Pringle et al. 1999). A 
pc-scale jet closely aligns in projection on the sky with the 
rotation axis (Herrnstein et al. 1997). 

Modjaz et al. (2005) present an analysis of polarimetric 
observations at 22 GHz of the water vapor masers in NGC 
4258 obtained with the VLA and the GET. They do not de- 
tect any circular polarization in the spectrum indicative of 
Zeeman-induced splitting of the maser lines of water, and ob- 
tained only an upper limit on the magnetic field strengths. 
They obtained the 1 — a upper limit value of the toroidal 
component of the magnetic field at a radius of 0.2 pc the 
value of 90 mC and determined a 1 — cr upper limit of 30 mG 
on the radial component at a radius of 0.14 pc. They also find 
from their observations of magnetic field limits that the geo- 
metrically thin disk model and the jet-disk model arc better 
candidates for accounting for the extremely low-luminosity of 
NGC 4258. 

More recently, Reynolds ct al. (2008) have shown from 
analysis of SUZAKU and XMM-Newton observational data 
that the observed iron lines originate in the surface layers of 
an warped accretion disk at the distance 10'^ — lO^i?^ from 
the black hole. In contrast to the majority of Seifert 2 galax- 
ies, there was no indication of a Compton- thick obscuring 
torus. The weak iron line and the lack of a reflection point to 
circumnuclear environment that is remarkably clean of cold 
gas. They note that such a circumnuclear environment is only 
found in two AGNs - NGC 4258 and M81 that contrast to the 



majority of Seifert 2 galaxies. 

This picture coincides with one by Herrnstein et al.(2005) 
pointed out earlier. According to them the observing intrin- 
sic absorption in the X-ray spectrum can arise in the outer 
layers of the warped geometrically-thin accretion disk at the 
distance ~ 29 pc from the black hole, where the molecular- 
to-atomic transition occurs. 

This picture allowed us to use our simple model of arising 
of optical polarized radiation without taking the warp contri- 
bution into account. Of course, the power-law dependence of 
the magnetic field is an important assumption in our deriva- 
tion. This assumption is now commonly accepted (see Pariev 
et al. 2003). 

The detected of polarized continuum and line emission from 
the nucleus of NGC 4258 was by Wilkes et al. (1995). After 
that, Barth et al. (1999) obtained spcctropolarimetric obser- 
vations of the NGC 4258 nucleus at the Keck II telescope. 
The observations were obtained on 1997 April 10 UT at the 
Keck II telescope with the LRIS spectropolarimeter. The re- 
sults of these observations are presented in the Table 1 of the 
paper by Barth et al. (1999). For the continuum polarization 
they obtained the following results: 

p(AA4000 4800 i) = 0.38 ± 0.03%, x = 12° ± 2° 

p{XX 5100 ^ 6100 A) = 0.35 ± 0.01%, x = 7° ± 1° 

p{XX 7500 8500 A) = 0.29 ± 0.02%, x = 8° ± 2° (16) 

We see that the polarization is weakly increasing to the short 
wavelength range, but the value of a position angle is practi- 
cally constant. The position angles x are the angles between 
wave electric field oscillations and the surface of the accretion 
disk. 

3.1 Estimates of magnetic field in the model 
of nonturbulent accretion disk 

For the inclination angle i = 83°, the expected polarization 
should have the value p(0,/i) = 6.9% (/i = cosi = 0.122). 
From Eqs. (9) and (10) for degree of polarization we find that 
possible parameters a = S\\ii and b = 6± \fY^^]? are near the 
value 20. The position angle x for this value of a is equal 
to 43.6°, which is far from observing values (see Eq. (21)). 
From the discussion in Sect. 2.1 we know that for high values 
of parameters a and 6, the possibility of small position angles 
exists if a ~ 6. 

The exact formulae (8) give us the values a = 122 and 
b = 122.9, which correspond to observed values of polarization 
degree 0.38% and x = 12° at A = OAAjim. These values 
correspond to 6\\ = 1000 (or B|| = 6400 G ) and S± = 124 
(or B± = 800 G). The analogous values of parameters for 
the polarization degree 0.35% and x = 7° at A = 0.56/zm 
are ^|| = 762.3 (or = 3037 G) and 6i_ = 95.5 (or Bj_ = 
380.5 G). For the case A = O.Sfim {p = 0.29%, x = 8°), the 
formulae (8) give = 2401.4 G and Bi_ = 298.3 G. We see 
that in all cases the normal magnetic field B|| is greater than 
B± . Because the effective temperature T,, decreases with the 
increase in the distance from the inner radius of an accretion 
disk, we conclude that magnetic field also decreases with the 
growing distance. 
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What seems unsatisfactory in these results is their sensi- 
tivity to small variations in parameters a and b. If these 
parameters change their values (A (6 — a) ~ 1) slightly the 
solution is impossible. For magnetic fields it means that the 
values have not to change its values greater than 10 G. This 
is very improbable for real situations in accretion disks. For 
this reason we have to seek a more satisfactory model where 
this sensitivity does not occur. Such a model really exists. It 
takes into account that the magnetic field can be turbulent. 

3.2 Estimates of magnetic field in the model 
of turbulent accretion disk 

According to Silant'ev (2005, 2007) the chaotic component 
B' of the magnetic field (B = Bq + B'), where Bq is a regular 
part of the magnetic fic;ld, gives rise to additional extinction 
of the intensity of linearly polarized waves (parameters Q and 
U) due to small scale chaotic Faraday rotations. The Gaus- 
sian distribution of turbulent velocities was assumed. Math- 
ematically, this effect is analogous to the known problem of 
diffusion of scalar impurity in a stochastic velocity field (see, 
for example, van Kampen 1981). In our case, the Faraday 
rotation angle replaces the role of impurity. The main part 
of the effective cross-section at, = arC corresponding to this 
additional extinction, takes a very simple form: 



C = fB n{6'^) ~ 0.64/s riA*(Mm)(B'2(G))/3. (17) 

Here, n is mean Thomson optical radius of turbulent eddies, 

/b ~ 1 is a constant characterizing a particular form of two- 
point turbulent velocity correlations (for estimations we take 
Jb = 1), and 6' and B' are fluctuating parts of Faraday de- 
polarization parameter 5 and magnetic field B, respectively. 
Clearly, the additional extinction should be proportional to 
the mean square of fluctuations of physical parameter 6. Be- 
cause 5 X^, the dimensionless parameter C ^ X*. 

The asymptotic formulae taking this effect into account 
have the same form as formulae (3-5) with the substitution 
(1 — fc/i) ^ (1 + C). (Remember that we consider the con- 
servative atmosphere with g = A; = 0.) In particular, the 
formulae (8) acquire the form 



- / d<P 
71" Jo 



7!" Jo 



(O) = (l + C7)Q(0,M)x 

(1 + C)2 + a^ + b^ cos2 $ 
[(l-f C)2 a2 + 62 cos2 $]2 - {2ab cos $)2 

{U)^aQ{0,fj.)x (18) 

(1 + C)2 + - 62 C0S2 $ 

[(1 + C)2 +a'^ + tP cos2 $]2 - (2a6 cos i>)2 ' 

Here parameters a and b are presented by formulae (7), 
where and B± denote the regular (mean) values of cor- 
responding magnetic fields, i. e. i?o|| and Bo±. The numer- 
ical calculations show that the relative polarization degree 
p(B, n)/p(0, /i) is also the symmetric function of parameters 
a and 6, as for the nonturbulent atmosphere. 

The parameter C is fairly large for our case, C ~ 10 -j- 20. 
In this case the term (2a6cos$)2 in denominators of formulae 



(18) can be neglected and we obtain the following analytical 
expressions: 

?,2 

tan 2x 



(l + C)(l + d)' (l + C)2+a2' 



MO,m) 



[(l + C)2+a2]y(rTd) 



'(1 + C)2 + 



(l + rf)2 



(19) 



These expressions are fairly exact up to a, 6 < C. If pa- 
rameters a or 6 are equal to zero, formulae (19) are exact. 
In particular, for the cases of pure normal regular magnetic 
field Bo± = (6 = 0), and pure perpendicular regular field 
Boll = (a = 0), they transform to the following exact for- 
mulae: 



p(B,n) 



p(o,m) 



tan 2x ■ 



1-hC" 



(20) 



p(B, n) 



p{0,H) 



X 



0. 



V(l + C)2+5i(l-/,2) 

In the limit case of nonturbulent magnetic field (C = 0), 
they coincide with the formulae (9) and (10). If the chaotic 
magnetic field B' is fairly high (Ti\^(pm){B''^) > 4i?o||.oj_) 
sec Silant'ev 2007), then the depolarization of radiation is 
due to additional extinction (17): p{B',n) = p{0,fi)/C ~ 
1/A'*(i3'2). If we suppose that {B'"^) = Const, then p{X)turb ~ 

The existence of new parameter C describing the level of 
magnetic field fluctuations makes the estimation of mean val- 
ues of _B[| and B± difficult. In our case this problem is sim- 
pler because the level of magnetic fiuctuations (parameter C) 
changes slowly with the variations in parameters a and 6. In- 
deed, for the pure normal magnetic field (6 = 0) we found 
a = 7.37, C = 15.5 for A = 0.44/xm; a = 4.77, C = 18.13 
for A = 0.56/Ltrn, and a = 6.56, C = 21.88 at A = 0.8/Ltrn. In 
the limiting case of high values of polarization degree (o = 6), 
these values are: a = 6 = 8, C=14;a = 6 = 5, C = 17.5, 
and o = 6 = 7, C = 21, respectively; i.e., parameter C is 
practically the same for these two cases. 

For this reason we calculate the a and 6 parameters us- 
ing the mean values for parameter C; i.e., C = 14.75,17.81 
and 21.44, respectively. This gives the values a = 8, 6 = 6 
(B|| = 423 G, B± = 39 G) for A = 0.44/xm; for A = 0.56A<m - 
a = 5, 6 = 3.5 (B|| = 163 G, B_i_ = 14G); and for A = O.S/xm 

a = 7, 6 = 5.5 (5|, = 112 G,B^ = 10.8 G). These values 
of magnetic fields are lowere than those for the case of non- 
turbulent accretion disk. But important is that they were 
derived without restriction a ~ 6. 

Let us estimate the level of magnetic fiuctuations taking the 
mean values for parameter C. We also assume that fs = ^ 
and Ti = 0.1. We do not know the real distribution of turbu- 
lent eddies in a turbulent accretion disk. The general picture 
of turbulence consists of cascade of eddies with different di- 
mensions. The small eddies with ti ^ 1 includes a value 
propto a parameter G. The large eddies do not allow us to 
describe the considered effect in the range of radiative transfer 
equations. It seems that our value n ~ 0.1 is fairly natural 
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for describing the turbulent effcc;ts in magnetized plasma. If 
we increase n to k-times, then the level of magnetic fluctu- 
ations (-B'^) decreases to k-times. The real value of parame- 
ter Ti can therefore be estimated by independent estimation 
of magnetic fluctuations. After that we obtain the values 
B' « 136 G, 92 G, and 49 G for the mean square root values 
of magnetic fluctuations at places where the thermal radiation 
has a maximum for A = 0.44/im, 0.56/im, and 0.8/xm, respec- 
tively. These values of fluctuations are equal to 31%, 56%, 
and 44% of the mean magnetic fields for mentioned wave- 
lengths. 

3.3 Estimates of the magnetic field at the 
horizon of the black hole in NGC 4258 

Wo now proceed to the estimation of the magnetic field 
strength in NGC 4258 using the maser polarimetric data. The 
magnetic field structures in accretion disks are difficult to ob- 
serve and remain poorly known. If the disk is penetrated by 
a dipole field of the central object or by a global field of the 
surrounding interstellar medium, there may be a net vertical 
flux. Sano et al. (2004) consider the models of an accretion 
disk with a uniform magnetic fleld. The stress forces in ac- 
cretion disks may be proportional to (Hawley et al. 1995) 
or to (Tmncr ct al. 2003). 

Zhang & Dai (2008) have studied the effect of a global 
magnetic fleld on viscously rotating and vertically integrated 
accretion disks around compact objects using a self-similar 
treatment. They show that the strong magnetic fleld in the 
vertical direction prevents the disk from being accreted and 
decreases the effect of the gas pressure. 

On the other hand, Konigl (1989) and Cao (1997) underline 
that the inclination of the fleld lines at the surface of the disk 
plays a crucial role in the magnetically driven outflow. They 
show that, for the nonrelativistic case, a centrifugally driven 
outflow of matter from the disk is only possible if the poloidal 
component of the niM,(i,n(^tic field makes an angle less than a 
critical 60° with the disk surface. 

Now let us estimate the value of R\ for the model of a 
standard accretion disk using the observational data of NGC 
4258. The spatial structure of standard accretion disk have 
been calculated by Poindexter et al. (2008). The size scaling 
is determined by Eq. (14). The basic physical parameters of 
the central nucleus of NGC 4258 are Mbh = 3.9 • IO'^Mq and 
?7 = Lboi/LEdd = 10-2-^'^ = 0.0054 (Satyapal et al. 2005). 
We use these estimations below in numerical calculations. 

The estimates of the scaleradius Rx from Eq. (14) give the 
following results: 



Rx = 1.38 • 10^^ ( — 



1/3 



cm, Xrest = 0.44 fim, 



0.1 



1/3 



Rx = 1.9 • 10 I — ) cm, Xrest = 0.56 /Ltm, 



0.1 



1/3 



i?A = 3.06 • lO^"" 1^— j cm, Xrest =0.8 ^im. (21) 

For £ = 0.1, which corresponds to u/Mbh — 0.6 in the model 
of Novikov & Thorne (1973), the ratios of Rx/Rh are equal 
to 



Rx{Q.U)/Rh = 13; Rx{QM)/Rh = 18; 

i?A(0.80)/i?H 29. (22) 

The formulae (21) for arbitrary wavelength A acquire the 
form: 



Rx = 1.94-10^ 



^A^ 

.1/3- 



(23) 



These estimates show that Rx < Rab, if the viscosity param- 
eter a > 0.1 (see Eq. (23)). Here Rab is the boundary radius 
between two zones of the standard accretion disk (Shakura & 
Sunyaev 1973): (a) a radiation pressure-dominated zone and 
(b) a gas pressure-dominated zone. In both zones the opac- 
ity comes from Thomson scattering. The boundary between 
these zones is given by Shakura & Sunyaev (1973) and Pariev 
& Colgate (2007): 



«(5^)"'(l^)"'(?)"' 

This means that the characteristic spatial radius of an ac- 
cretion disk Rx, corresponding to observed wavelength A, lies 
for NGC 4258 in the zone (a). For the frequently used value 
a = 0.01 (see Pariev & Colgate 2007), we have Rab = 23.2i?i/ 
at e = 0.1; i.e., the inequality Rx < Rab takes place only for 
A = 0.44/xm and 0.56/im. 

We use the polarimetric data by Modjaz et al. (2005). 
These data allow a 1 — cr upper limit of Bmas ^ 30mG on the 
radial component of the disk magnetic field at the radius of 
0.14 pc. Using the power-law radial dependence of magnetic 
field (11), we obtain the following expression: 



Br, 



B{Rx) 



( Rx 



5/4 



= 0.018 1^ 



5/12 



(25) 



where R^as = 0.14 pc, B{Rx) = 425 G for A = O.Unm. 
This expression can be used for a crude estimate of parameter 
e. In our case we obtain e « 0.03. In this case, Rab = 
58.2i?ff(a/0.01)2/2i. This means that inequality Rx < Rab 
occurs even for value a = 0.01, mentioned above. The values 
Rx for e = 0.03 are higher than values (22) to 1.5 times, i. e. 
they also lesser than Rab, if we take a = 0.01. 

Finally, we can estimate the magnetic field strength Bh 
at the horizon radius of the black hole in NGC 4258 using 
the data for Rx/Rh presented in expressions (22), and for 
our values B{Rx) for the turbulent accretion disk model (see 
Sect. 3.2). Taking in Eq. (11) r = Rx, we obtain the expres- 
sion: 



B 



H 



5/4 



(G). 



(26) 



As a result, we obtain for A = 0.44/im, 0.56fim, and A = 
O.Sum the following values Bh = 1.05 • lO'^G, 6.06 • lO^G and 
7.57 ■ lO^G, respectively (at e = 0.1, which corresponds to 
spin (Kerr) parameter u/Mbh — 0.6 in the model of Novikov 
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Table 1: The value of Bh [G] for various data of Kerr param- 
eter and radiative efBciency. 



Source 



£ = 0.081 



: = 0.998 
£ = 0.32 



Mbi 



= -0.9 
£ = 0.039 



NGC 4258 41; 3.5 6.6 • 10^ 



4- 10^ 



9.6 • 10^ 



Table 2: The value of linear polarization from accretion disks. 



Source log 



Mb 



log 



Lsdd 



p 



0.5 



Mbh 

e = 0.081 



P [%] 
0.99S 



Mbh 

e = 0.32 



BLRG 


8.64 


-1.5 


0.06 


0.057 


NLRG 


8.14 


-2.92 


0.09 


0.067 


SSRQ 


9.29 


-0.90 


0.05 


0.057 


FSRQ 


9.01 


0.44 


0.94 


0.042 


Seyfert 1 


7.23 


-0.59 


1.40 


0.036 



& Thornc 1973). At e = 0.03, these values are to 1.66 times 
higher. This case approximately corresponds to a /Mbh ~ 
-0.95 (sec Krolik 2007). 

Our estimates are slightly different as a result of different 
error intervals of polarimetric data (see Eq. (16)). Besides, 
some uncertainty exists in the choice of the level of fluctua- 
tions (parameter C). It seems for A = 0.56/im that this uncer- 
tainty is less than for other wavelengths. For this reason the 
estimate Bh = 6.06 • lO'^G seems to be preferable. As an ex- 
ample, we present the estimates for this effective wavelength 
and other various parameters s and a /Mbh in Table 1. 

In the estimations, presented above we take n = 5/4 in 
the basic Eq. (11). How do we change the estimations for 
other values of n? The calculations give the following values 
of Bh at the event horizon of the supermassive black hole in 
NGC 4258: = 3 • 10^ G at n = 1, and Bh = 5.3 • 10^ G at 
n = 2. These results indicate that the magnetic field strength 
of SMBH in NGC 4258 at the event horizon should be at the 
level « 10^ - 10^ G. 

It should be noted that our estimations do not use the 
values of viscosity parameter a. It was only shown that, for 
a > 0.1, the radius Rx lies inside zone (a). Of course, our 
estimations depend on parameters e and the power-law index 
n, which, in principle, are to be found in a detailed model of 
an accretion disk. 

The data in Pariev et al. (2003) are given for r > 2.5rg 
{rg = 2GMbh/(? is the gravitational radius of black hole), 
i. e. slightly beyond the horizon radius Rh (beyond 2.2Rh 
for £ = 0.1, and beyond 5Rh for e = 0.03). It seems that 
they used the simplified theory, which do not "work" near 
the horizon. The calculations in Pariev et al. (2003) corre- 
spond to magnetic dominated regime; i.e., magnetic energy 
is greater than radiative thermal energy. Because R\ lies in- 
side the zone (a) (radiation dominated zone), we derive that 



plasma parameter /3 = Pgas/Pmagn < 1 in model of Pariev 
et al. (2003). Here Pgas and Pmagn are gas and magnetic 
pressures, respectively. 



4 Magnetic coupling process in 
AGN and QSO: testing by contin- 
uum polarization 

Li (2002), Wang et al. (2002, 2003), Zhang et al. (2005) 
studied the magnetic coupling process (MC) as an effective 
mechanism for transferring energy and angular momentum 
from a rotating black hole to its surrounding accretion disk. 
This process can be considered as one of the variants of 
the Blanford-Znajek (BZ) process (Blanford & Znajek 1977; 
Blanford & Payne 1982). It is assumed that the disk is sta- 
ble, perfectly conducting, thin, and Keplerian. The magnetic 
field is assumed to be constant on the black hole horizon and 
to vary as a power law with the radius of the accretion disk. 

Since the magnetic field on the horizon Bh is brought and 
held by its surrounding magnetized matter of a disk, the some 
relation must exist between the magnetic field strength and 
accretion disk and, finally, the bolometric luminosity of AGN 
(see Ma et al. 2007) . This relation takes a form 



Bh{G) = 



k^/\2L,oi/ecY/^ 



R 



h 



(27) 



Here Li^oi = £Mc^, Rh is the horizon radius (see Eq.(12)), M 
is the accretion rate, e the radiative efficiency (calculated, for 

example, by Novikov & Thornc 1973; Krolik 2007; Shapiro 
2007, see also Sect. 2.2). Coefficient k is the inverse plasma 
parameter [3 = Pgas/Pmagn, where Pgas and Pmagn are gas 
and magnetic pressures, respectively. We shall consider the 
case /3 = 1 {k = 1) in future, i.e. the equipartition case. 
Then Eq. (27) is transformed into 



Bh{G) 



Me 
M. 



BH 



1/2 



(?) 



1/2 



6.2 • 10^ 



(a/M, 



(28) 



BH) 



where q/Mbh is the Kerr parameter, 77 = Lhoi/LEdd and 
LEdd = 1.3 • 1 0^^ (Mbh /M0) (erg s~^) is the Eddington lumi- 
nosity. 

We next calculate the expected polarization value of radia- 
tion in a number of specific AGNs, taking the effect of Faraday 
depolarization into account, as considered above. In mech- 
anisms of magnetic coupling, we mainly consider the case 
where magnetic field is supposed directed along the normal 
to the accretion disk (see, for example, Wang et al. 2002, Ma 
et al. 2007). Our estimates of the magnetic field in NGC 4258 
(see Sect. 3.2) show that the vertical magnetic field is much 
larger than the perpendicular one. For these reasons, we also 
consider only B|| fields as first approximation, so we have to 
determine B\\ = B{Rx) from Eq. (25). The values Rh and 
Rx are presented in Eqs.(12) and (14), and Bh - in Eq. (28). 
As a result, we obtain the following formula for dimensionless 
parameter ^m: 



0.8A2(/xm)B||(G) ~ 474A^/3(/xm) 



1/12 



M. 



BH 



Mr. 



-1/12 



1/4 



(29) 
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The estimations of polarization degree p{X) can be obtained 
from Eq. (9). 

A systematic analysis of a large sample of AGN available in 
the BeppoSAX pubhc archive was performed by Grandi et al. 
(2006) . Their sample includes AGN of various types. Narrow 
line radio galaxies (NLRG), broad line radio galaxies (BLRG), 
steep spectrum radio quasars (SSRQ) and flat spectrum radio 
quasars (FSRQ) (see Table 6 from their paper, where the 
values Mbh /Mq and ij = L^oi / Ledd are presented. 

The results of our calculations of the dimensionless depo- 
larization parameter 5\\ are 

BLRG : Sn ~ 79; 60; NGRG : Sn ~ 67; 57, 



SSRQ : ~ 79; 60; FSRQ : ~ 108; 82, 



Seyfert 1 : (Jii ~ 125; 95, , 



(30) 



where the first numbers correspond to £ = 0.081 and 
a /Mbh = 0.5, and the second ones correspond to s = 0.32, 
a/ Mbh = 0.998. 

The estimations of possible p{^) for the inclination angle 
of the accretion disk i = 60° (p = cosi,p(0,/i) = 2.25%) 
are presented in Table 2. Because all values of » 1, we 
can use the approximate formula p{iJ,) ~ p(0, /u)/5||/i. This 
formula is valid up to i = 87°, so for i = 85° (p(0,/i) ~ 
7.8%, fj, = 0.087) the polarization degree is to 20 times 
higher than values presented in Table 2 (p(/i = 0.087) = 
(p(0, 0.087)/p(0, 0.5)) (0.5/0.087) ~ 20). For other inclination 
angles i, the calculations are analogous. We stress that po- 
larization observations of sources, presented in Table 2, do 
not exist up to now, so we present only possible values of 
polarization degrees. This procedure, presented below, can 
be used to estimate the source parameters if the polarization 
data arc available. 

The limiting case a/Msn = corresponds to e = 0.057; 
i.e., the predicted polarizations differ only slightly from the 
presented values for u/Mbh = 0.5, so for the source NGC 
4258 instead of 6.6-10'^ G, corresponding to u/Mbh = 0.5 (see 
Table 1), we obtain 7.6 • 10^ G for limiting case u/Mbh = 0. 

5 The wavelength dependence of po- 
larization of AGN and accretion 
disk models 

Polarization in AGNs can be intrinsic or extrinsic. Light scat- 
tering by a nonspherical distribution of electrons near the 
central engine of AGN is a basic intrinsic polarization mech- 
anism; namely, an accretion disk is a typical example of a 
nonspherical distribution. Scattering by magnetically aligned 
dust grains in the interstellar medium of galaxies is the typ- 
ical example of an extrinsic situation. The very important 
feature characterizing the polarized radiation from a magne- 
tized accretion disk is the wavelength dependence of polar- 
ization degree that is very different from that of Thomson's 
scattering. 



For a strong magnetic field strength, when 6\\^± ^ 1 the 
simple asymptotic formulas follow from Eqs.(9) and (10): 



f(A) 



p(o,m) ^ p(o.m) 
p(o>m) ^ p(o,m) 



(31) 



For turbulent atmospheres, with the existence of chaotic mag- 
netic field component B' , the corresponding asymptotic for- 
mulae follow from Eqs. (20). If the parameter of turbulence 
C ^ 6 ^ 1, then p{X) l/A'^, if the level of magnetic fluc- 
tuations {B'"^) is constant (see Sect. 3.2). Below we mainly 
consider the limiting spectra (5||,_l 3> 1) only for nonturbulent 
atmospheres. For turbulent atmospheres with a high value of 
C, the corresponding spectra are to be divided to A^. 

The wavelength dependencies of the radiation flux and its 
polarization essentially stem from the radial distribution of 
the temperature in an accretion disk. For a standard accre- 
tion disk (Shakura & Sunyaev, 1973), radial dependence of 
the temperature takes the form: ^ r~^/*. To get the inte- 
grated spectrum from the disk, we add up all of the Planck 
curves from each radius. If ~ then the radiation flux 
(see Pringle & Rees 1972, Shakura & Syunyaev 1973, Gaskell 
2008): 



(32) 



For a standard accretion disk (s = 3/4), the flux (30) takes 
the known form Fi, ~ u^^^. 

Substitution of Tg ^ into formula (13) leads to rcla- 

1 /s 

tion Rx ~ Ae between the characteristic scale Rx and the 
effective wavelength of polarimetric observations. According 
to expression (11) we obtain B\\^^{Rx) ^ A~"/^. As a re- 
sult, from Eqs. (31) and (11), we obtain the next wavelength 
dependence for the case of a strong Faraday depolarization: 



P(A) 



ff(0,M) 
5||,±A2 



A"/ 



s-2 



(33) 



As a result, the degree of observed polarization p{X) depends 
on both power-law indices n and s. This comes from the 
depolarization effect of Faraday rotation of the polarization 
plane. In the absence of this effect, the polarization degree, 
independent of radiation flux Fx, depends only on the inclina- 
tion angle of the accretion disk. Of course, the usual Stokes 
parameters I,Q, and U give the information on the radial 
distribution of the effective temperature (e.g. Eq. (32)). 

For distribution of magnetic field in a standard accretion 
disk (s = 3/4) with n ^ 3/2 {B\\ j_ ~ i?"^/^), the polarization 
does not depend on the wavelength. For a standard accretion 
disk with n — 5/4, mostly elaborated by Pariev et al. (2003), 
the wavelength dependence of the polarization is quite weak 
p{X) ~ A-i/3. In this case the wavelength dependence of the 
Stokes polarized flux takes the form p{v)Fj, ~ i/^/^. In the 
general case the formulae (32) and (33) give rise to relation 



(34) 
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The presented formulas allow us to test the various models 
of an accretion disk using the data of the wavelength depen- 
dence of polarization of AGN and quasars. 

For the turbulent accretion disk with a high value of pa- 
rameter C, formula (34) transforms to the expression 

p{u)F^ ^u"^-"^. (35) 

The position angle values Xj used in our formulae, denote the 
angles between the accretion disk plane and the observed di- 
rection of electric field wave oscillations. If it is impossible 
to estimate the disk inclination, wc cannot obtain the the- 
oretical X from the observed position angle. The observed 
spectra of polarization degree po&s(A) and the differences of 
position angle Xobsi^i) — Xobsi^j) are to coincide with the 
corresponding theoretical values p{X) and the differences in 
position angles x{X)~x{\) for different wavelengths Xi and 
Xj. These conditions make the estimation of magnetic field 
and unknown disk's inclination angle more precise than for 
the single wavelength observations. 

5.1 Discussion of some observational data 

Webb et al. (1993) and Impey et al. (1995) present the data 
of measurement UBVRI polarizations of a sample of AGNs 
and QSOs. The position angle appears to be wavelength- 
independent, suggesting that the polarization in a given ob- 
ject originates in a single physical process. In many cases 
the percentage of polarization increases with frequency. Au- 
thors have compared the polarized fluxes with the predictions 
of competing models of polarization in AGNs: synchrotron 
emission, scattering from electrons or different types of dust 
grains, and electron scattering in an accretion flow. In nine 
sources from this sample, the polarization seems to be the 
result of dust grain scattering. A number of these sources 
(NGC 4151, Mrk 509, NGC 5548, Mrk 290) has best charac- 
teristics of model due to electron scattering in an accretion 
disk or torus. It is interesting that, for these four sources, the 
slopes of F„ - ly'f are close (7 ~ -1.37,-0.9,-0.86,-0.9, 
respectively; see Webb et al. (1993). According to for- 
mula (32), we obtain the corresponding values of parameter 
s = 0.46, 0.47, 0.52, 0.51: i.e., they are near .s = 0.5. 

The general shape of the polarization spectrum with pa- 
rameter "q" was determined by a least-squares fit propor- 
tional to p{v)F,^ = Av'^ . It appears that the slope of the 
wavelength dependence of polarized flux q varies widely, be- 
tween —2 and with a typical uncertainty of 0.3. 

What values of parameter n corresponds to limit values of 
q, if we accept s = 1/2 using formulas (34) and (35)? For 
nonturbulent accretion disk (Eq. (34)) we obtain 



This formula gives n = 1.5 for q ~ —2, n = 1 for q = — 1, 
n = 0.5 for g = 0, and n = for g = 1. Remember (see Pariev 
et al. 2003) that accretion disk with decreasing gas density 
at large distances r corresponds to n > 1. It seems the self- 
consistent models of Pariev et al. (2003) allow the values 
of n ~ 1 — 2. The high values with n > 2 correspond to 
a very sharp decrease the magnetic field inside the accretion 



disk. That is why we investigate the cases 1 < n < 2 as 
physically acceptable according to the results of Pariev et al. 
(2003). From Eq. (32) we see that px ~ A^^""!) at s = 1/2. 
This means that the degree of polarization that diminishes 
with the increasing A cannot be explained inside the physical 
limitations, assumed above (n > 1 corresponds to increasing 
polarization degree with A). 

For this reason we investigate what gives the very high 
magnetic turbulence. For this case, formula (32) transforms 
to expression p{X) ^ A"/*~^. It leads to p{X) ~ A^<^"~^^ at 
s = 1/2. As a result, our physically acceptable restriction 
1 < n < 2 allows us, in principle, to explain the decreasing 
degree of polarization. Taking the definition p{v)Fi, ~ v'^ into 
account, we derive from Eq. (34) the following relation 

n={7 -q)s- 2. (37) 

For s = 0.5 this expression gives n = 1 for q' = 1, and n = 1.5 

for q = 0. For q = —2 this formula gives n = 5/2, which 
characterizes a very sharp decrease in magnetic field inside 
the accretion disk. 

As an example, we now consider the continuum radiation of 
source NGC4151. According to Webb et al. (1993), we have 
7 ~ —1.37, and from Schmidt & Miller (1980) the parameter 
q is equal to (7 ~ —0.33 (see Fig. 3 for the value p{v)Fi, in 
the interval A = (0.3 - 0.94)/xm). From Eq. (31) we obtain 
s = 2/(3 — 7), which gives in this case s ~ 0.46. Substitu- 
tion of the values q = —0.33 and s = 0.46 to Eq. (37) gives 
n ~ 1.37. In nonturbulent case we have n ~ 0.45, which cor- 
responds to the increase in gas density with growing distance 
from the nucleus (see Pariev et al. 2003). In our case, when 
wc know the observed spectra ~ v'^, and p{u)Fi, ~ v'^, 
both Eqs. (34) and (35) give rise to the same expression 
p{X) - AT-9 ~ A-104. Thompson et al. (1979, see Fig. 3) 
present the polarization degree p(A)% in the same interval 
A = (0.3 — l)/zm. The data have a rather wide distribution. 
Using our dependence, we can approximate the polarization 
degree by the formula p{X) ~ —0.5 -|- 0.7 ■ X^^ '^'^{jjLm), which 
fairly satisfactorily describes the mean value of the presented 
data. 

6 Conclusions 

We have presented the method for estimating the magnetic 
field strength at the event horizon of a supermassive black 
hole through the polarization of accreting disk emission. The 
polarized radiation arises from the scattering of emission light 
by electrons in a magnetized accretion disk. Due to Faraday 
rotation of the polarization plane, the resulting polarization 
degree differs essentially from the classical Thomson case, be- 
cause the wavelength dependence of the polarization degree 
appears. This feature means that the magnetic field strength 
at the event horizon of a black hole can be estimated from 
polarimetric observations in the optical range. 

For estimating observed polarization, we use the az- 
imuthally averaged asymptotic expressions of the Stokes pa- 
rameters for outgoing radiation, assuming that the accretion 
disk is optically thick and that the Milne problem takes place. 
Using these formulae, we discuss the wavelength dependence 
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of the observed spectra of polarization degree. We also con- 
sider turbulent accretion disks when the magnetic field pos- 
sess both regular component and chaotic one. Since the polar- 
ization spectrum of scattered radiation strongly depends on 
the accretion disk model, our results can be used to construct 
a realistic physical model of the AGN environment. 

The estimates of the magnetic field strength of supermas- 
sive black holes in NGC 4258 are presented. We also found 
for the source NGC 4151 that the observed polarization de- 
gree spectrum can be satisfactorily explained by our mecha- 
nism with the high level of turbulent magnetic field. In this 
case the power-law regular magnetic field is B • 



[11] Dolginov, A. Z., Gnedin, Yu. N., Silan'ev, N. A., 1995, 

Propagation and polarization of radiation in cosmic me- 
dia, 

(Gordon & Breach, N.Y.) 



-1.37 



-0.46 



and 

In the cases where the observed polarization is 
the result of various mechanisms, such as light scattering in 
an accretion disk and a jet, our method can be considered as 
important part of the problem as a whole. 
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